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Abstract 

One of the challenging aspects of electroproduction at high-energy is the un- 
derstanding of the transition from real photons to virtual photons in the GeV^ 
region. We study inclusive electroproduction on the proton at small xb using a 
nonperturbative dipole-proton cross section calculated from the gauge invariant 
gluon field correlators as input. By quark-hadron duality, we construct a photon 
light cone wave function which links the "hadronic" behavior at small to the 
"perturbative" behavior at large . It contains quark masses which implement 
the transition from constituent quarks at low to current quarks at high Q^. 
Our calculation gives a good description of the structure function at fixed energy 
for < 10 GeV^. Indications for a chiral transition may already have been seen 
in the photon-proton cross section. 

1 Introduction 

Since the late 60's, multi-GeV electron and muon collisions with protons have been inten- 
sively used in order to get information on the proton structure. For photon virtualities Q 
far below the Z-mass, the interaction is mediated by a virtual photon and the inclusive 
photon-proton cross section can be described by means of the proton structure functions. 
At large Q^, the leading unpolarized structure function, F2(xb, Q^), is to leading-log ac- 
curacy the well-known linear combination of partonic distributions, qf{xB, Q"^), weighted 
by the square of the parton electromagnetic charge expressed in units of the proton 
charge: 

F2{xb,Q^) ='^ejxBqf{xB,Q^). 
f 

An illustrative partonic description emerges when the process is envisaged in a frame 
where the proton has a large momentum p (formally |p| oo) and in a particular 
gauge. At small Q^, the above decomposition and the partonic interpretation of the 
process get spoiled by power corrections. 

There is, yet, an alternative to the infinite momentum frame description of the colli- 
sion which is a description in the center of mass frame. In this frame, the photon acquires 
a structure and we have to deal with the interaction of two structured objects. Although it 



may look as if we had not gained anything by changing our point of view, the operation is 
interesting if we focus on the high-energy fixed kinematical domain of the process. In 
this regime, the bulk of the photon-proton interaction ressembles that of hadron-hadron, 
i.e., diffractive scattering. In the Regge approach which is applicable in the kinematical 
region envisaged here, this is understood as being due to the universality of the Pomeron. 
Implicitely, this assertion assumes that the photon has developed an internal structure 
due to its coupling to strongly interacting quark fields and that this structure gives the 
main contribution at small xb to its interaction as compared to the direct contribu- 
tion from its "bare component" . Therefore a great deal of insight can be gained from a 
common understanding of both photon-hadron and hadron-hadron collisions. 

In Ref. the application to diffractive scatterings of hadrons of the model of the 
stochastic vacuum has been carried out. Recently, in Ref. 0, the same approach was used 
to describe diffractive leptoproduction of vector mesons in the range = 2-10 GeV^, 
thus starting to implement the program just mentioned. Our aim in the present paper is 
twofold: we want to pursue the comparison by considering the total photon-proton cross 
section and we want to extend the phenomenology to 0. 

In Ref. , the interaction amplitude for the exclusive vector meson photoproduction 
off a proton has been written as 

-^rv^,{z,r)j'^\z,r,t). (1) 

il)v and V'-y are, respectively, the vector meson and virtual photon light cone wave func- 
tions. If the final vector meson wave function is replaced by the virtual photon wave func- 
tion, one gets the forward Compton amplitude, that is iso"*"^. The quantity J^^\z,r,t) 
represents the Pomeron exchange amplitude for scattering of a qq dipole of size r off the 
proton target, where an average over the dipole orientations has been carried out. The 
light cone fraction carried by the quark in the photon is denoted by z. The Mandelstam 
variables for the process are s = W"^ and t. The photon is characterized by its virtuality 
Q"^ and polarization. 

The quantity J^^-* has been derived in the model of the stochastic vacuum |0 following 
the method of Ref. 0. In these references, the few parameters which fix the magnitude 
and shape of J^"-* have been adjusted to fit the phenomenology of proton-proton elastic 
cross section. We use the same parametrization here and focus on the photon structure. 

Let us discuss the photon wave function that enters in Eq. (||). There are two standard 
schemes. The first one is to expand the photon wave function in a hadronic basis 
The wave function for a transversely polarized photon reads 

i^%=E^H^,MT) + Rest, (2) 

p,L0,4> y ~r V 

where we have explicitly written the low-lying vector meson contribution. The symbol 
Rest stands for a sum over residual 1~ states, like higher radial and orbital excitations 
and nonresonant multiparticle states. According to the vector meson dominance (VMD) 
hj^othesis, the low-lying vector meson states, p, u, 0, dominate the photon wave function 
at small making Eq. (^) a useful expansion in this regime. The wave function for a 
longitudinally polarized photon is obtained by changing efyMy — > efvQ, ipv(T) ''Pv(l) 
and the Rest term accordingly. 



2 



To assess the relative importance of the residual contribution Rest at large Q^, we 
examine first the experimental behavior of the inclusive cross section and the ratio = 
<7l/<7t in the range = 1-10 GeV^ (cf. left part of Table 0). At large Q^, the success 
of the parton model with spin 1/2 quarks tells us that the structure function F2{x) oc 
Q^(j{Q'^) scales, i.e., the total cross section decreases as and that it is dominated by 
transverse photon scattering. In VMD the contributions of p, u, alone would lead to 
a 7*-p total cross section ctt oc Q~'^ and ctl oc Q~^, i.e., a dominating longitudinal cross 
section. It means that in the VMD description the transverse inclusive cross section must 
be built from the residual term Rest in the photon wave function. Let us next consider 
vector meson production (cf. right part of Table |1]). In the range = 1-10 GeV^ the 
experimental cross section has a behavior and is predominantly longitudinal. This 
dependence is overshooted by the longitudinally dominated VMD cross section 0(Q~^). 
It thus turns out that the Rest term in the photon wave function Eq. (0) is also needed 
to cancel the p, u, contribution in order to provide the right behavior for the 
elastic vector meson production. This cancellation mechanism has to be implemented in 
a consistent way to get a unique description of the photon in high energy scattering. 







exp 
VMD 
free qq 


< 0.2 > 1 and rising 

g-2 g2 Q-2 Q2 
g-2 const g-6 g2 



Table 1: Large behavior of the cross section a = ai + (Tt and of the longitudinal to 
transverse ratio R = ctl/ctt for both total photon-proton cross section and elastic vector meson 
production. The first line shows the experimental results in the range = 1 — lOGeV^. The 
second line displays the asymptotic behavior expected from VMD. For comparison, the third 
line gives the scaling behavior determined with a free qq wave function (up to logarithms). 



A possible solution is to use a quark-gluon basis. At values of large enough, 
the hadronic component of the photon is indeed a free qq pair, if)^^'^ k, ip^q^ so that the 
quark-gluon basis is more efficient than the hadron basis and the cancellation mechanism 
occurs automatically. It seems phenomenologically possible to envisage a kind of hybrid 
description where the photon state can be viewed as a superposition of a few low-lying 
resonances plus a free gg-state. We stress, however, that identifying the Rest in Eq. (H) 
with the free qq wave function ipqq is not sufficient because, on the one hand, it would 
not reproduce the phenomenology given in Table |I|, and, on the other hand, it leads to a 
double counting of some hadronic configurations. In order to avoid these problems, one 
has schematically to modify Eq. (0) in such a way that, for small values of g^, it agrees 
with the vector-dominance like form 

and, for large g^, it approaches the perturbative photon wave function. 

This is the problem for which we propose a solution in the present paper. A central 
role in our parametrization of the photon wave function will be played by the effective 
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quark mass. With increasing resolution of the photon the hght quarks experience a sort 
of chiral transition with constituent masses at low becoming current quark masses 
for >1 GeV^. The outline of the paper is as follows. In Sec. |^ we demonstrate that the 
two dimensional harmonic oscillator can be used to model our light cone parametrization 
of the photon wave function. In Sec. ^the approximate form of the photon wave function 
is given at low Q^. Section ^ deals with the calculation of the inclusive virtual photon 
cross section for large energy and small x^. In Sec. ^ we discuss corrections from finite 
energy. 

2 The two dimensional harmonic oscillator as a model 
for the photon wave function 

In light cone perturbation theory, the photon wave function is given by the light cone 
energy denominator and spin matrix elements. Leaving aside the spin complexities, we 
have the approximate form: 

r (P'\^ e^^^ 1 

Here mj is the current quark mass of the quark and antiquark with flavor /. The trans- 
verse extension of the wave function is given by r± ~ e~'^, where e = ^ z{l — z)Q'^ + m^. 
For small values of e, however, the confining gluonic forces and/or the spontaneous chi- 
ral symmetry breaking will intervene and limit the transverse extent of the photon wave 
function. At large energy, far away from the target, the perturbative wave function is 
certainly valid, but at finite small Q^, there is enough time for the photon to dress up 
like a bound state. 

In quantum mechanics, the two-dimensional harmonic oscillator is a very reasonable 
testing ground for the behavior of the photon wave function in transversal space, since 
the harmonic oscillator has two essential features in common with the behavior of our 
qq dipole in QCD: on the one hand, large transverse distances are prohibited, because 
of the harmonic potential (confinement), and, on the other hand, the potential vanishes 
at the origin which corresponds in our problem to the fact that, for short times and 
small relative transverse distances of the quark and antiquark, the dynamics is entirely 
governed by the kinetic energy in the hamiltonian (asymptotic freedom). The Green 
function of the two dimensional harmonic oscillator 

G(r,0,M)= E , ^-^^y^f^ (4) 

n={ni,n2) V ' ^ ' ' 

shows the analogy to the photon wave function. The wave function '?/'n(0) stands for 
the hard process of qq production and V'n(r) gives the transversal extension. The short 
time restriction can be included by looking at the dynamics for large negative values 
oi E = — M, where large M corresponds to the deep Euclidean region of QCD. The 
harmonic oscillator Green function approaches for large negative values oi E = —M the 
free two-dimensional Green function in quantum mechanics 

72 1, ^ikr 



d'k e m 



GU^, 0, M) = J = -K„iV^ (5) 
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When we put 2mM = —2mE = z{l — z)Q'^ + m'j, we see directly the similarity to 
the perturbative photon wave function. Since many exact results are available for the 
harmonic oscillator, we will be able to check several manipulations on the transverse 
wavef unctions. In the end, we will not apply the results of the harmonic oscillator di- 
rectly to QCD, but we extract the essential features from the non-relativistic model and 
transpose them into relativistic quantum field theory. 

To begin with, let us illustrate the point made in the introduction on the cancellation 
mechanism. In diffractive vector meson production, we need the matrix element of of 
the qq state with the vector meson state. In the harmonic oscillator case, this second 
moment, 

00 



J rf^r r2G(r,0,M)^oo(r), (6) 

can be computed exactly for the full Green function. Using the spectral decomposition of 
the full Green function, Eq. (Wj, and the one-dimensionnal harmonic oscillator property 



(n|x2|0) 

one gets 

r') = — 
1 00 iom 



nO 



72^^02(0) V^oo(O) 
3tu + M w + M 



In the last equation, the contribution of the second excited state and the ground state 
are shown separately. With '?/'oo(0) = y^cum/vr and V'02(0) = —^uj-mjl-K^ they add up to 



^ Im \l 7Tm{3uj + M){uj + M)' 
The free Green function gives for the same moment 

As expected, the large M behavior of the full Green function moment is exactly repro- 
duced by the free Green function. The important lesson of this simple exercise is the 
demonstration that the large M behavior follows from a delicate cancellation between 
ground and excited state contributions. Note that the contribution of one single state is 
~ 1/M and would overshoot the full result at large values of M. This is to be compared 
with the discussion we had in the introduction. 

Although the harmonic oscillator Green function, G(r, 0, t), is known analytically, we 
know of no such representation for the Fourier transformed G{r, 0, M = —E); therefore, 
we have obtained an "exact" expression by performing the sum in Eq. with the 
first 500 terms using the known wave functions of the harmonic oscillator. A simple 
calculation shows, that even for moderate values of M, say 5uj, one needs more than 
twenty intermediate terms in the representation Eq. (|^) in order to obtain a better 
accuracy than the one from the free Green function. In terms of vector dominance, 
this means that we need many intermediate vector mesons in order to get an adequate 
description for moderate values of the photon virtuality. As we will show, a much more 
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efficient procedure is to shift the argument 2mM in the free Green function by a M 
dependent value Sq] it turns out that this method gives, even for M = 0, a very decent 
approximation to the full Green function. 

Since the exact Green functions is available for the harmonic oscillator, the shift 
parameter can be calculated by comparing the modified free Green function with the 
exact one. This procedure cannot be, however, carried over to the lightcone wave function 
of the photon. We consider, therefore, the "two-point" function, n(M), and its derivatives 

Jn Jn 

nW(M) := (-l)"^n(M) := (-1)«^G(0, 0, M), (7) 

instead of the "three-point" functions G{r,0, M). For convergence, we need at least one 
derivative. The two-point functions have been extensively studied in QCD, especially 
with sum rule techniques @. From this, we know that in an asymptotically free theory 
the ansatz "one resonance plus perturbative continuum" is a very good phenomenological 
representation for the two-point function in the Euclidean region. Our method is then to 
make for the two-point function, 11*^"^) (M), the model "one resonance plus perturbative 
continuum" and to use for the three-point function an approximate form which can be 
parametrized easily and adjusted in such a way that the two-point function obtained 
from it agrees with the model two-point function. 

For the derivatives of n(M), the ansatz "one resonance plus perturbative continuum" 
reads, to lowest order, 

rW..^^_ ^!|V^oo(0)P , m (n - 1)! 



(cj + M)"+i 27r(st/2m + M)"^ 



where ipoo is the ground state wave function and St the continuum threshold above which 
we use the perturbative Green function. Duality states that the integral from to St 
over the imaginary part of the free (i.e., perturbative) two-point function accounts for tt 
times the residue at the resonance pole: 

rfA;2lmnfree(A;2) = vr|V^oo(0)p. 



Thereby, we obtain the continuum threshold 

St = 4:7i\iJoo{0)\^ = ium. 

As a simple approximate Green function, we consider the free Green function with a 
shifted argument 

TTl i 

Ga(r,0,M,so) = -iro(V2mM + sor). (9) 
The derivatives of the correponding approximate two-point function have the form 

ni")(M, So) = i-ir^^GM 0, M, So) = — , . (10) 

(Notice that in the left hand side the relation involves the partial derivative and not 
the total derivative.) By equating the approximation Eq. ([To|) to the phenomenological 
function Eq. (||) 

nW(M,so) = n(^(M), (11) 
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we can determine the only free parameter of the approximation, namely so(M), and 
obtain 

So/2m = {uj + M)("+^)/"(2cu + M)[2nuj{2uj + M)" + {u + M)"+^]-^/" - M. (12) 

The exact form of the shift depends on the number of differentiations assumed. Yet, at 
M = 0, so/2m is around uj/2 and decreases to become a small correction to M, i.e., less 
than 5%, for M larger than 3uj. 

In Fig. |I|, we display the approximated Green function obtained with a n = 3 shift 
(long dashes) and the full Green function (full line), for M/u = 0, 0.5, 1 and 4. For the 
first two values, we also show for comparison a two-resonance approximation to Eq. (H). 
For the last two values, the short dashes represent the free Green function; for M/u = 4 
the approximated Green function can hardly be distinguished from the free one. 




0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 




0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 



Figure 1: Green functions (in units of m) of an harmonic oscillator as function of r (in 
units of l/y/um) for different values of M/uj. Solid line: exact Green function G{t, 0, M), 
long dashes: our approximation G'a(r, 0, M, sq), i.e., the shifted free Green function Eq. (^ 
with the shift sq from Eq. (p^; short dashes: approximation with two resonances, dots: 
free Green function. 

We estimate the quality of our approximation by forming the moments 

(r'^) := j dh r''G\r,0,M), 

which can be evaluated easily, both for the full and the approximate Green function. Since 
in electroproduction, for high virtualities, the gg-proton cross section behaves approxi- 
mately like and, in the hadronic region, like r^'^ 0, we pay special attention to the 



7 



moment (r^). In Table we give the relative differences of the exact and approximated 
moments 

A2 := ((r^)-(r\)/(r^), 
Ao := ((r°) - (r°),)/(r°), 

for different values of M and numbers of differentiations n. As can be expected, the lower 
n- value, n = 1, yields a good approximation for the zeroth moment, whereas the second 
moment is well reproduced with n = 3. The maximal error for (r^) turns out to be 11% 
(at M=0). 



M/uj 


n - 


= 1 


n 


= 3 


n - 


= 5 




A2 


Ao 


A2 


Ao 


A2 


Ao 





1.09 


0.01 


0.11 


-0.26 


-0.16 


-0.35 


0.5 


0.51 


0.02 


0.05 


-0.14 


-0.13 


-0.23 


1 


0.28 


0.02 


0.02 


-0.09 


-0.11 


-0.15 



Table 2: Differences of the exact and approximated moments: A2 := {{t"^) — (^^)a)/('"^) 
and Ao := ((r") - (rO),)/(rO) 

We also compare the overlap of the exact Green function and the ground state wave 
function, i.e., (?'^)oo defined in Eq. (|^), with the respective overlap including the approx- 
imated Green function 

('■')ooa = / ^'G'a(r,0,M)^oo(r). 

This is shown in Fig. |^ and Table ^. One sees that the shifted free Green function gives 
an good estimate of the exact matrix element, the relative error being at most 10%. 



M 


approx 


exact 


error 





0.400 


0.376 


0.06 


1 


0.155 


0.141 


0.10 


2 


0.0819 


0.0752 


0.09 


3 


0.0504 


0.0470 


0.07 


4 


0.0342 


0.0322 


0.06 


5 


0.0247 


0.0235 


0.05 


6 


0.0186 


0.0179 


0.04 


7 


0.0146 


0.0141 


0.03 


8 


0.0117 


0.0114 


0.03 


9 


0.0096 


0.0094 


0.03 


10 


0.0081 


0.0079 


0.02 



Table 3: Comparison of exact and approximated overlap with ground state (in units 
m = a; = 1). 
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0.001 0.01 0.1 1 10. 

M 

Figure 2: Overlap Eq. of the Green function of a two-dimensional harmonic oscillator 
with the ground state; solid line: exact result, dashes: calculated with our approximated 
Green function from Eq. @ (in units m = = 1). 

In Table ^ we compare the shifts sq obtained from Eq. ([12|) with n =2, 3 and 4 
differentiations, with the ones needed to get exact agreement for the second moment. To 
reproduce the second moment, the displacement sq is quite different from the value of sq 
which one obtains from Eq. (|ll|). Our method underestimates the shifts considerably for 
large M, but the overall errors on the matrix elements remain small as shown in Table |^ 
and Fig. 0. Indeed we already noticed that for M ^ u the full Green function Eq. @ 
is well approximated by the free one Eq. (H), i.e., sq can be safely set equal to in this 
region. 



M/iu 


n = 2 


n = 3 


n = 4 


exact 





0.485 


0.547 


0.593 


0.585 


1 


0.279 


0.344 


0.397 


0.455 


2 


0.179 


0.234 


0.283 


0.385 


3 


0.123 


0.168 


0.211 


0.325 


4 


0.090 


0.127 


0.162 


0.27 


5 


0.068 


0.099 


0.129 


0.255 



Table 4: so/2m determined by the sum rule method for n =2, 3, 4 subtractions and 
adjusted to give exact overlap. 
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3 Approximate photon wave function extended to 
low values of 

We now want to apply the approximation methods developped for the harmonic oscillator 
to the photon wave function. For this purpose, we consider first the polarization tensor 
for the vector current = ijj'y^i{j of a quark of mass m 

W{q,m^) = J d^xe"^''{0\T[J^{x)r{0)]\0) 

At large g^, the imaginary part of the polarization function n(g^, m^) is obtained to 
lowest order in perturbation theory from the free quark-antiquark propagation. One has 



2 2 + 2m^ Am? 

Imn(g ,m ) = ^^^f ' (13) 



The polarization function itself is only determined up to a subtraction constant, but its 
derivatives 

nW(Q2 = -q^m") := -^n{q\m'), (14) 
can be written for n > 1 by dispersion relations: 

n(n).^2 _ T , Imn(.,m^) 

Due to asymptotic freedom, the polarization tensor has a very good representation 
in the Euclidean region consisting of the ground state vector meson with mass my and 
residue Fy and perturbative qq continuum calculated with current quark masses: 



V is the lowest lying vector meson in the flavor channel considered, i.e., p,uj = {uu =F 
dd) I \/2 and = ss. Fy is the decay constant of this vector meson deflned through 

Fymv£^(g,A) = (0|J^(0)|\/(g,A)). 

The continuum threshold St can be related to the decay constant Fy by local duality: 

Fl = - r dslmU{s,ml). 

TC J Am? 



f 



We copy now the procedure from the discussion of the harmonic oscillator. The photon 
wave function plays the role of the three point function which we want to approximate. 
As in Eq. (P), we take as approximate Green function the free Green function but shift 
the variable Q^. The structure of the perturbative photon wave function is of the form 
(see Eq. (|)): 

ijj^ oc Kq {\Jz{l - z)Q'^ + m|r) . 
10 



A dependent shift thus corresponds to a replacement of the current mass m/ by 
an effective mass wicfr(Q^)- Here some improvement is still possible, since in reality the 
virtuality appears in combination with the light cone momentum fraction through the 
term z{l — z). We leave this difficulty aside for the moment and shall return to it later. 
From Eqs. (|T^-(|T3p, we form the derivatives of the approximate polarization function 

^t\Q\ml^) ■■= ^(^^^,(Q^m,\) (16) 



127r2 (9(-g2) 



Next, in complete analogy to Eq. (0), we determine the effective mass in such a way 
that the shifted two point function, ni"\ is equal to the model two point function H^^: 



n\F^ n\ Imn(s,m^) 

(g2 + ^2^)n+i + 74 ^(7Tg2 

The method gives the effective quark mass as a function of on the left hand side 
from the purely hadronic parameters Fy and my on the right hand side. In Fig. |^, we 
display the resulting effective quark masses mesiQ"^), putting n = 2, for massless current 
quarks and strange quarks with a current mass = 150 MeV. As for the harmonic 
oscillator, the result depends on the number of differentiation assumed. For rriq = 0, the 
effective mass starts at = at 220 MeV for n = 2 and 245 MeV for n = 3 which are 
typical constituent quark mass values. The effective mass drops to at Qo ~ 1.5-2 m^. 
Beyond this value of Q^, it formally becomes imaginary. We have seen in the harmonic 
oscillator that the method underestimates the shift at large virtuality, but the errors 
introduced are less than 5%. We therefore put the effective mass equal to zero above Qq. 
To be specific, we use the simple linear parametrization 

meff(Q') = .22 (1. - Q^Ql), in GeV, for < Ql = 1.05 GeV^, 
meff(Q') = 0, for > q2_ ^^g) 

To support our point to set the effective light quark mass to zero for large Q^, we display 
in Fig. I the second derivative of the model Green function together with the free one. We 
see that the two agree at large values where our effective mass formally would become 
imaginary. A more refined procedure would be to make a smooth connection between the 
small behavior obtained with the present method and the behavior obtained around 
1 GeV^ using Operator Product Expansion [^]. We shall not dwell on this possibility in 
the following. For the strange quark, we refer to the second part of Fig. ^, which shows 
that the corresponding starting value for the constituent strange mass is 310 MeV. It 
reaches its asymptotic value in the range = 1.5-2 GeV^. A simple parametrization 
for the strange quark mass is 

rrisesiQ^) = .15 + .16 (1. - Q^Ql), in GeV, for < Ql = 1.6 GeV^, 
m,efr(Q') = .15, for > q2_ ^^g) 
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Figure 3: Effective quarlc masses which reproduce the model polarization function, see 
Eq. (p^. (a) The behavior for light quarks and (b) strange quarks. 



As mentioned before, in the photon wave function, appears together with the 
factor z(l — z). We re-analyzed the polarization function 11^'^) including an effective 
quark mass as a function of Q\q = Az{l — z) . This modifies the treatment that leads 
to Eq. ([I6|) along the following lines. The imaginary part of the polarization function 
n(g^,m^) of the vector two point function Tl^^i, can be written as 

Im 11( g , m ) = / az H z[l — z)q — m 

12tx Jo q 

Upon integration over z, it yields to the familiar result given in Eq. (p!3|). We then obtain 
the derivatives by dispersion relations: 

Setting m = niesiQes) the right hand side of this equation gives a new expression for 
n^"). We numerically find that the phenomenological two point function can be repro- 
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Figure 4: Second derivative of the polarization function; dashes: model function nph(Q^), 
i.e., one resonance plus continuum; solid: lowest order perturbative expression with zero 
mass. Note for > 1 GeV^ the perturbative expression becomes indistinguishable from 
the perturbative one. 

duced with an accuracy better than 10% (see Fig. ^ with parametrizations similar to 
that given above, provided we change the value of Q^. We get for light quarks 

mefr(Qe'ff) = -22 (1. - Qes/Qo ), in GeV, for Ql^ < Q'^ = 0.69 GeV^, 
meff(Q,\) = 0, for Ql^ > Q'^ (21) 

The lower scale Q'q is due to the difference between the average in z occuring in Eq. ( pOf ) 
and the value at the average z = 1/2. For the strange quark, the result is 

m.efr(Qeff) = -15 + .16 (1. - Qls/Q'i^), in GeV, for Ql^ < Q'^ = 1.16 GeV^, 
m.efr(Qe\) = -15, for Q^s > Q'l (22) 



4 Inclusive photon-proton cross section 

Let us now consider the forward Compton amplitude, i.e., the amplitude in Eq. (0), with 
the replacement ipy ^ ip*, aX t = 

M{j*p^-1*p) =is j ^^^^\ip^{z,r)\'^ Jp{z,r), 

where we use the notation Jp{z,r) = J^^\z,r,t = 0) for short. We employ the same 
dipole-proton cross section as in our previous work on exclusive vector meson produc- 
tion 0, which is based on the evaluation of gluon field strength correlators between 
Wilson areas mapped out by the color neutral dipole and the proton. The absolute size 
and z, r dependences of the cross section Jp{z, r) are determined by the gluon condensate 
(g'^FF) = 2.49 GeV^, the correlation length a = 0.346 fm and the transverse radius of 
the proton Rj_p = 0.52 fm, together with the form of the correlators assumed there. 
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Figure 5: Second derivative of the polarization function for vector current of light 
and strange quarks. Solid: model function nph(Q^), dashes: approximate expression 
na(Q^,meff) with the quark mass depending on Ql^ = Az{l — z)Q'^ . Note that for the 
quark mass depending on alone we get exact agreement by construction. 

We compute the square of the photon wave function using the expression given in 
Ref. P]. This leads to the longitudinal and transverse amplitudes 



pi roo 

Ml/t = isJ2(^fL/T = isj^*^} rdrlL/T{z,r) 

J y «/ 

with the integrands 



Uz\l-z)''Q''Ko{eryjp{z,r), 
{[z^ + (1 - zY]e''K,{ery + m}KoierY} Jpiz, r). 



47]-2 
47]-2 



(23) 

(24) 
(25) 



The extension parameter of the photon is e"^ = z{l — z)Q'^ + 'm'j. It depends on the quark 
flavor through the quark mass and thus each flavor contributes in a different way to the 
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above sums. Kq, Ki are modified Bessel functions. Tliey arise from the perturbative liglit 
cone wave function of the photon when one takes into account all the spin complexities 
ignored in Sec. 0. For a general dipole-proton cross section, Jp, these expressions are 
identical to those given in Refs. [H, |^. Our dipole-proton cross section Jp{z,r) calculated 
from gluon-gluon correlators describes the scattering on a proton of a loop with transverse 
size r and infinite extension along the light cone. It depends only very weakly on the 
longitudinal momentum fraction z and shows for values of r smaller than approximately 
2a ^ 0.7 fm the dipole behavior 



Jp{z, r) ^ Cr with C 



4.3 



(26) 



For very large values of r it increases linearly and around 1 fm it goes approximately like 



a. 5 



One can get a rough idea on the dependence of the cross sections if one assumes 
the dipole behavior given in Eq. (^). Then one obtains analytical expressions for the 
cross sections for each flavor channel, setting u = Am? /Q"^: 
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From these expressions, we deduce easily the behavior for large and small Q^: 
1) large Q"^: 



CTfL 
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ln-^-- + — In-^ 



(27) 
(28) 



2) small Q 



2. 



(TfL 



.2 2N, C Q2 

D 



l + 0{Q^/m^ 



l + OiQ^m' 



(29) 
(30) 



For large values of and longitudinal photons the small dipole sizes are dominant 
and the use of a purely perturbative photon wave function is justified. We see indeed that 
the limit of the longitudinal cross section for large values of is unproblematic. For 
transversal photons however small values of the longitudinal momentum fractions z and 
(1 — 2;) become more important and large dipole sizes may play an important role even if 
the value of is large. A behavior of Jp{z, r) thus leads to a logarithmic divergence 
in the quark mass, as can be seen from Eq. (0). For large values of r such a behavior is 
however unrealistic, and any reduced increase of the form r^"*^ with e > leads to a finite 
cross section even for a perturbative photon wave function and a quark mass equal to 
zero. This feature stresses the importance of a realistic model both for the long distance 



15 



part of the photon wave function and the dipole cross sections. A good qq wave function 
and a reahstic dipole-proton cross section become even more relevant for low values 
of Q^. Our model for the hadronic scattering part is inherently nonperturbative. For 
the photon we absorbe the nonperturbative effects in a virtuality dependent constituent 
mass. We stress that all input to these models is taken from sources outside the realm 
of electroproduction. 

4.1 Photoproduction 

Let us compare our computed cross section with data, starting with photoproduction. 
In the following we consider data at a center of mass energy W = 2Q GeV. We choose 
this value because it is the one where the model parameters are adjusted to fit the corre- 
sponding proton-proton elastic scattering data. At this center of mass energy, the photon- 
proton total cross section is 118 /ib. The Pomeron part of the Donnachie-Landshoff fit |0 
gives cTpom = 110 /ib, the remaining part being attributed to other Regge trajectories. 
The model of the stochastic vacuum accounts for the Pomeron part of the cross section. 
The Reggeon contribution to strange (and heavier) quark interaction is much suppressed 
(Zweig suppression) and can be neglected. The charm quark contribution is measured 
independently and is rather small, 1 /ib. In the present study, we focus on u, d and s 
which give the bulk of the cross section. According to Donnachie and Landshoff, the 
strange quark contribution to the total cross section is 8.3 /ib. 

Our theoretical results for the light, u and d, and strange quark contributions are: 

au+d = 84 /ib, 
as = 9.4 /ib. 

Some comments are called for. To illustrate our purpose, it is again useful to stick to 
the approximate amplitudes derived by assuming the short distance behavior Eq. (pG]). 
In Eq. (|2^) the magnitude of the cross section, given by C ~ 4.3, is determined by 
the properties of the QCD vacuum, namely the gluon condensate and the correlation 
length of gluon field strength correlators, and the proton radius. As we have seen in the 
beginning of this section, the transverse extension of the photon is bounded by the value 
of the quark mass, r < l/rries, therefore the dependence of the cross section is certainly 
correct if the quark mass is large enough, otherwise it gives a first approximation. The 
limit of small was given in Eq. (|29|) : 

This means that the photoproduction cross section depends crucially on the value of the 
constituent quark mass. The reproduction of the cross section within 15 % is encouraging 
especially taking into account the fact that we have determined the effective quark mass, 
Eq. ([18D , without any recourse to electroproduction phenomenology. If we attributed the 
remaining 15 % difference between the experimental value of the cross section and our 
theoretical one to the value of rries, we would get a 8 % decrease of the effective quark 
mass value. Due to other sources of uncertainty, this refinement does not make sense 
here. We also notice that, in the approximate cross section, 40% come from the second 
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term in Eq. (|25|) which is proportional to the quark mass squared. Within the present 
determination of rries, this supports the interpretation of the modification of the photon 
extension parameter as being due to the generation of the effective quark mass rather 
than being just a shift in the argument of Ki in Eq. (^Sf). For the strange quark, the 
comparison with the extracted cross section is correct within 10%. As an indication, we 
note that our amphtudes with a current strange quark mass of 150 MeV would produce 
a much too big cross section of 31 /xb. 

4.2 Electroproduction 

We now consider virtual photon, ^ 0, scattering off a proton. We form the structure 
functions 

= —^i^T + ^l) (32) 
Fl{Q') = —^OL. (33) 

Since the light and strange quarks contribute in a different way, we first calculate the 
corresponding quantities separately. Special attentions will be paid to the dependence 
of the structure functions at fixed W = 2^ GeV, which corresponds to the energy where 
we determined our input dipole-proton cross section. We want to investigate the question 
whether the chiral transition from the constituent quark to the parton can be seen in the 
inclusive electron scattering data. In our theoretical calculation the effective quark mass 
mefT(Q^) evolves with the photon virtuality . A priori it is not clear whether the photon 
virtuality or the combination of light cone momenta and Q^, namely = 4z(l — z^Q^ ^ 
should be used in the running of the quark mass for low virtualities. In the second case the 
integration over the light cone momentum fraction z may wash out the chiral transition 
effect. 

In Fig. H, we show the theoretical results for F^'^'^ixB = Q'^/W'^,Q'^) and -F|(^s = 
Q^/W^, Q^) at a fixed energy W = 20 GeV. At this energy and in the range consid- 
ered, the Reggeon contribution to F2 is less than 5% |^. 

The calculation in Fig. |^ with fixed light quark mass tuq shows a dependence similar 
to the logarithmic dependence ln((5^/mo) expected for the dipole-proton cross sec- 
tion. This agreement is also quantitatively good. The sliding quark mass m = rricsiQ'^) 
produces a steeper variation of F2 for 0.1 GeV^ < < 1 GeV^, in agreement with 
the variation of the mass in the logarithm. Above = 1.05 GeV^ where its value has 
gone to zero, the finiteness of F2 is connected with the long distance behavior of the 
dipole-nucleon cross section. In this region of the z averaged integrand in Eq. (^Sj) 
has a maximum in r in the region r = 0.2-0.8 fm and it extends to very large values of 
r, being only damped because of the change of the quadratic behavior of Jp = Cr"^ at 
small distances to a power 1.5-1.8. The maxima of the integrands for the transverse and 
longitudinal cross sections in r are not so different but the profile is much faster decreas- 
ing in the longitudinal case. In the case of a dependence of nics on = 4:z{l — z)Q^ 
the resulting structure function F2 in Fig. |^ interpolates smoothly the case of mesiQ"^) 
between the minimal and maximal values of Q^. It lies about 10-20% above the curve 
with constant mo for > 0.2 GeV^. The strange quark structure function reaches at 
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Figure 6: (a) Theoretical calculation of F2 as a function of . The long-dash curve is 
the light quark contribution with the effective quark mass given in Eq. ([I8|). The solid 
curve is with the effective quark mass given in Eq. (^). The short-dash curve is with a 
fixed constituent quark mass mo = 0.22 GeV. The strange quark contribution with the 
effective mass given by Eq. ([T9|) is shown as a dashed dotted curve. The dotted curve is 
with the effective strange mass of Eq. (p2D. (b) Same study for the longitudinal structure 
function Fx,. 
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Figure 7: = ai/cr-r as a function of Q^. The solid curve is our expectation with z 
dependent quark masses for hght and strange quarks, Eq. (^Tj) and Eq. (^), respectively. 
The dashed curve is the result with the effective quark masses given in Eq. ([T8|) and 
Eq. Data are from NMC 

the maximum the asymptotic rate of 20% of the light quark structure function. The 
difference between the contribution given by the effective mass of Eq. (|T9|) and the one 
of Eq. (^21) is qualitatively the same as in the light quark case. 

The longitudinal scattering is suppressed compared with the transverse scattering in 
agreement with the estimates. Here the difference between the z dependent effective 
quark mass and the only dependent quark mass is very small, since the longitudinal 
photon wave function suppresses the endpoints in the 2;-integration. Conversely the dif- 
ference with the fixed uiq mass reaches about 30% at = 1 GeV^ and is more visible 
than in F2. 

At the photoproduction point the ratio R = (TL/aT vanishes, it then increases until 
Q2 = 1 GeV^. For the quark mass depending on the virtuality its behavior is more 
flat than when the quark mass depends on Qg^. A computation with a fixed mass, mo, 
and the short distance behavior Jp oc leads to a ratio similar in shape with a maximum 
in around 40 mg. This result was already obtained in the classic paper of Bjorken, 
Kogut and Soper on the electroproduction of lepton pairs in a slowly varying external 
field 0. Note that in this case a constant quark mass mo = 0.22 GeV would give a 
maximum in i? = (Tl/ctt at ~ 2 GeV^. The precise behavior of the ratio R = ai/ ut 
gives another signal for the chiral transition in the experimental deep inelastic scattering 
data. We show in Fig. ^ the ratio we get in our computation, combining light and strange 
contributions. We compare this ratio to NMC results. Unfortunatly experiments do not 
reach the small transition region. 
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Figure 8: Contribution of u, d and s to F2 as a function of Q^. Curves are as in Fig. |^. 
Squares are NMC results H and diamonds are E665 [ITU . 



We next combine the light and strange contribution to form the structure function 
F2 = F^^'^ + F| as a function of and for fixed energy W = GeV. In Fig. |, 
we compare our results with data from both E665 and NMC at energies 18.5 GeV< 
W < 21.5 GeV. The transition region is correctly described and the scheme with the z 
dependent effective mass seems to be preferred. At large region, our result overshoot 
data by 10-20%. As we now show finite energy corrections may reduce cross sections in 
this large region. 



5 Modification at finite energy 
5.1 Slow quarks 

In the forward Compton amplitude, we integrate over configurations with quark lightcone 
fractions varying from to 1. For the transverse cross section, there is a large contribution 
from aligned jet configurations, where one quark carries most of the momentum and the 
other one a minute fraction. At large but finite energy, however, the slow quark may not 
carry enough energy to generate a hadronic final state. Formally, the photon-proton total 
cross section is related to the imaginary part of the forward elastic amplitude via the 
optical theorem. In this amplitude we must sum only over the accessible channels, i.e., 
we have to take into account energy conservation at finite energies. To be more precise 
we have to care about how the energy is distributed in the physical color neutral final 
states. We require that the intermediate quark-antiquark and quark-diquark states (see 
Fig. 1^) both have an invariant mass bigger than a typical mesonic or baryonic state with 
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mass Mm = {Mp,MK*) or Mb = {Mn,Ma) respectively. 



{q-k + lf 



zi{l - Z2)W^ > Ml 
{l-z,)z2W''>Ml. 



(34) 
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Figure 9: The intermediate rearrangement for light quarks is determined by the momenta 
of the coUison partners as calculated in Eq. (|^) . 

For the present discussion, we restore the internal degrees of freedom of the proton, 
i.e., the cross section envisaged here is written as in Ref. with the full dependence on 
transverse and light cone coordinates of the quarks in both the photon and the proton. 



a. 



TP 



dZ2d?TC2 

Air 



|V'p(2;2,r2)|V(xi,Xi,X2,X2)e(zi,2;2, ly), 

(35) 

where 1 and 2 refer to the photon and nucleon sides, respectively. The proton is considered 
in a simple quark-diquark picture which has proved to be a good approximation in the 
framework of the stochastic vacuum model. In Eq. (|35[), we have inserted the threshold 
factor B which realizes the requirement of Eq. (|3^ ) : 



e{zi,Z2,W) = 9 Z,{1-Z2)W 



Mi 



9 



[1 - zi)z2W^ - M 



M 



(36) 



In the integrand of Eq. (p3D, the effect of the threshold factor is to generate a z = zi 
dependent phase space factor: 



<l>(z) 



dz2f{z2)9 



1 - 



Z2 



9 



Z2 - 



Ml 



'l-z)W^ 



In the integral, f{z2) represents the Z2 dependence as it results for the various term in 
Eq. (|35D. Indeed, since the Z2 dependence of the quantity J is rather weak, f{z2) turns 
out to be given essentially by the square of the proton wave function, for which we use 
f{z2) = 252z|(l — Z2Y- For large energies, W 3> Mm + Mb, ^{z) = 1 for intermediate z 
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and rapidly decreases to when z approaches the boundaries which at large enough W 
read 

Ml Mh 
— - < z < I — 

- - w^- 

Mp and = 20 GeV, the boundaries are approximately 0.0022 < 



Mp, Mm 



For Mb 
z < 0.9985. 

The threshold factor is therefore important in case where the endpoint contribution 
is sizeable. Usual hadron wave functions suppress this region in both hadron-hadron 
collisions and vector meson production, thus rendering this effect unimportant at large 
energy, s > 100 GeV^. In inclusive photon-hadron scattering, however, this endpoint 
region cannot be overlooked. The importance of the various region in z, in the full integral 
Eq. (|23|) , may be studied by varying the lower limit Z of the integration over quark light 
cone momenta in the photon wave function. Recall that the threshold condition for light 
quarks gives a lower limit z > 0.0022 for 1^ = 20 GeV as calculated above and that this 
lower limit goes like l/W"^. Since the integrand is symmetric in z, the upper integration 
limit can be restricted to z < 0.5. The resulting function It{Z) 



It{Z) 



0.5 roo 

dz / rdrlr(z,r) 
z Jo 



is shown in Fig. |TU| for the transverse cross section. 
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Figure 10: Interaction amplitude, IriZ)/ It{0), at = 3 GeV^ as a function of the cut 
off Z limiting the momentum fraction of quarks in the photon. The effective mass used 
in this computation is the one given in Eq. (|21|). 



In order to understand what happens in the endpoint region, it is instructive to study 
the behavior of the amplitude for a simplified Jp{z, r) behaving like C r" with n = 1 or 2. 
Our actual Jp is in some sense interpolating between those two choices, the power n = 2 
corresponding to the short distance dipole behavior already mentioned in Eq. (^61) . With 



22 



such a simple dependence, one can perform the r-integral analytically: 



^^^^^ = ^[.(l-.)g2 + m^]iW2' 

(rrlrT = r A^ + Vn)[z' + ^ ' , .oox 

y ^ [2(1 - 2)Q2 + ^2]n/2 [^(l-2)Q2 + ^2]l+n/2' ^^^^ 

where the overall constant is 

iV,C^,_, r^(l + n/2) 

For <^ Q^, one can focus on the first term in Eq. (pS]) . At Z ^ 0, the quantity 
Q'"'It{Z) behaves, for = 2, like 

Q%(Z) oc 0(1) - ln[Z + m^/Q\ 

and, for n = 1, like 



Q/t(^) oc 0(1)- J^ + mVQ 



In Fig. |I0], we see the transition from the logarithmic behavior given in case n = 2 to a 
linear behavior at very small Z. The limit Z = is however approached like in the case 
ra = 1: it does not have the dramatic dependence on exhibited by the n = 2 case. For 
< 4m^, the integrands become more or less fiat in z and the Z behavior is linear in 
the whole range. 

Because the photon size parameter, e"^ = z{l — z)Q^ + m^, is the important scale 
in X, the behavior changes in the neighborhood of Z = m? /Q"^. The region where the 
threshold suppression in photon-proton collisions is sizeable is therefore given by 



/g2 < M'^/w"' 



For M = Mp and W = 10 GeV, this shows that the effect becomes sizeable when 



mVM^ < xb- 

and W = 10 GeV, 

> 0.5 GeV^, for an effective quark mass ~ 0.1 GeV. If we were to consider a current 
quark mass below this value of Q^, the effect would show up at even a much smaller 
value of Q^. 



5.2 Threshold effect in the cross section 

We now compare with data our computed cross section modified by the threshold effects 
from slow quarks. For photoproduction, the change is negligible. For electroproduction, it 
is best to consider fixed and vary the energy W to see the threshold effect. In Fig. |1^ we 
show the variation of F2 for 9 GeV< W <25 GeV separately for = 1, 3 and 9 GeV^. 
(For convenience we shifted downward the first two sets by respectively 0.25 and 0.15.) 
We have checked that for these values of the inclusion of the Reggeon contribution 
given in Ref. |^ only modifies weakly the trend of our curves for W > 10 GeV. In each 
case we present results combining u, d and s contributions for both schemes: z dependent 
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Figure 11: Contribution of u, d and s to F2 as a function of W ^ for =1, 3 and 9 GeV^, 
from bottom to top (0.25 has been substracted to the first set and 0.15 to the second). 
The sohd hues represent our result for z dependent effective masses and the dashed hues 
correspond to dependent effective masses. Squares are NMC results and diamonds 
are E665 fig. 



quark masses, Eqs. (^) and (^2|), and z independent quark masses, Eqs. ( ]T8| ) and 
The effect is relatively stronger for the latter scheme since the endpoint contributions are 
clearly more important in this case. The difference between the two schemes decreases at 
small W where only intermediate z are taken into account. As anticipated the threshold 
effect is also more important at large and aX = 9 GeV^ the suppression of the 
cross section is typically 10% at IV = 20 GeV and reaches 30% aXW = IQ GeV. At this 
high the effective quark mass has practically gone to zero in both schemes and the 
large distances are cut off by the threshold condition. 

To complete our study we give in Fig. our result for F2{xb = Q^/W'^,Q'^) at 
W = 20 GeV including the threshold effect. The plot may be compared with Fig. |^ 
and we see that the qualitative aspect of the transition from low to high remains. In 
general the large range is sensitive to the threshold effect whereas the small range 
tests more the effective quark mass. 



6 Discussion and Conclusions 

We have computed the total photon proton cross section in a model of nonperturbative 
QCD. For values of the virtual photon mass > 2 GeV^ our input is the treatment 
of two Wilson loops in Minkowski space-time within a special model of nonperturbative 
QCD which approximates the infrared behavior by a Gaussian stochastic process deter- 



24 



F2 



0.1 - 



0.5 - 



0.4 - 



0.3 - 



0.2 - 




0.1 



10 (GeV2) 



Figure 12: Contribution of m, d and s to F2 as a function of Q^, for W =20 GeV. Curves 
and data are as in Fig. ^ 

mined by a non-local gauge invariant gluon field correlator. The latter one is essentially 
given by the local gluon condensate and the correlation length. In order to fix the size 
distribution of one loop a proton valence quark wave function has to be introduced. 
In principle all parameters of the model can be determined by sources other than high 
energy scattering, namely lattice gauge calculations and low energy phenomenology. In 
practice the errors in the parameters still necessitate some adjustment to high energy 
scattering data. In our case we have chosen the determination of Ref. based on the 
proton-proton total cross section and the logarithmic slope of the elastic pp cross section 
at zero momentum transfer. Although the parameters turn out to be rather stable there 
remains still a certain variation in a range of a few percent for the correlation length 
and 20-30% for the gluon condensate and the proton radius. The model gives energy 
independent cross sections, contrary to the slight s^'^^ energy dependence seen in low 
experiments. Since the parameters are related to the ISR energies, 20 GeV, we also 
compare the photon cross sections calculated here with experimental values around that 
energy. We have already used this procedure in order to calculate diffractive electropro- 
duction of vector mesons with very good results. 

In the present work we have extented the range down to = 0. In order to do 
that we have constructed nonperturbative photon wave functions essentially by intro- 
ducing a virtuality dependent constituent quark mass. We were encouraged to such a 
simple procedure by model investigations of harmonic oscillator Green functions. In our 
method we adjusted the value of the momentum dependent quark mass to reproduce 
the phenomenological two point function of the vector currents. The development of a 
quark mass at large distances has been seen also in calculations with the instanton liq- 
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uid model |TT|. We repeat that none of our input parameters was in any way related or 
adjusted to electroproduction phenomenology. 

Our approach differs in several ways from other investigations. For the treatment 
of the soft exchanges, Nikolaev and Zakharov have adopted a phenomenological two- 
gluon-exchange model [0. Their treatment also gives importance to the qq photon wave 
function and their dipole-proton cross section has a behavior at short distance and 
saturates in the r = 1-2 fm region. To suppress the contribution from large distance in 
the photon wave function they cut off the large size component using a smooth gaussian 
exp(— r^/i?^) with a confinement size parameter R,, k, 1.5 fm. In Ref. a current quark 
mass mu^d = 10 MeV was used, but in a later study of the same authors on the BFKL 
pomeron [jl2l a constant light quark mass 171^4 = 150 MeV was considered. As we have 
seen in our study such a value can help to limit the extent of the photon wave function, 
but at large there is no reason that the light quarks have masses different from the 
current quark masses. In the region of transverse distance probed, we have seen that 
the nonperturbative features of the gluon correlators are important and we think that a 
perturbative two-gluon exchange model cannot be trusted. 

Concerning the transition to small Q^, we have shown an approach different from 
vector meson dominance (VMD) frequently used in the low range. We claim that 
with increasing one would have to put a rapidly growing number of resonances into 
the VMD model which thereby becomes untractable. Our scheme of quark-hadron duality 
exploits the knowledge about the residue and the mass of the lowest vector meson state 
contributing to the vector current two-point function. 

Based on the transition between a VMD description of F2 at small and the par- 
tonic one at large Badelek and Kwiecihski have proposed to represent the proton 
structure function via dispersion relation |T3. The construction adopted here shares 



similarities with their approach although the latter makes no connexion to the notion of 
wave function which we need for the microscopic description of diffractive scatterings. 
In their approach the partonic contribution to F2 is extracted from structure function 
analysis at large and thus naturally fulfills perturbative QCD evolution at large vir- 
tualities where this evolution is experimentally observed. Such perturbative corrections 
are not implemented in our approach. 

Our results are encouraging. Without any adjustment they agree with experiments 
over the full range from to 20 GeV^ within 10-20%. For the high range we have 
shown how finite energy corrections may account for at least a part of the discrepancy 
without changing any of the model parameter and without affecting the good description 
of the transition region where the modification of the perturbative qq wave function is 
at work. It is rather difficult to decide with present data what is the best scheme in 
our approach, i.e., whether or not the effective mass depends on the quark light cone 
fraction. In order to put our result in a larger perspective, we show in Fig. |1^ the curve 
one deduces from Fig. [T^ with the z independent masses, Eqs. (|^) and (|19D , if we adjust 
the normalization by —10% (which would amount, e.g., to a 5% reduction of the gluon 
condensate). A rescaling of the curve in Fig. | by —15% shows the same pattern. We 
describe almost perfectly the dependence of the data in the whole range examined. 
This points toward the possibility that the chiral transition is already seen in present 
data, the kink in the data at = 1 GeV^ being related to the vanishing of the quark 
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Figure 13: Contribution of u, d and s to F2 as a function of Q^, for W = 2Q GeV. The 
curve is the result shown in Fig. [1^ and for the dependent quark masses of Eqs. (|18|) 
and (|19D (dashed curve in Fig. |12]) multiphed by 0.9. This rescahng would correspond to 
a 5% decrease of the gluon condensate. Data are as in previous figures. 

mass at that value of (chiral restoration). The photon-proton cross section could also 
be easily adjusted by a slight decrease of the constituent quark mass. 

Because of the fundamental importance of chiral symmetry breaking for hadron 
physics, the question of chiral symmetry restoration with increasing virtuality of 
the photon deserves to be studied in more detail. More systematic data for F2{xb,Q'^) 
and Fl{xb,Q'^) at fixed W and varying Q'^ would help in deciding whether a marked 
change occurs between the low momentum domain and the perturbative domain in in- 
elastic electron scattering. It may be that chiral symmetry restoration is seen more easily 
in electroproduction at a given virtuality Qq than in heavy ion collisions at a finite tem- 
perature 27lT = Qq. 
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